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ABSTRACT 

We evaluate the success of linear tidal-torque theory (TTT) in predicting galactic- 
halo spin using a cosmological TV-body simulation with thousands of well-resolved 
haloes. The proto-haloes are identified by tracing today's haloes back to the initial 
conditions. The TTT predictions for the proto-haloes match, on average, the spin 
amplitudes of today's virialized haloes if linear growth is assumed until ~ to/3, or 55- 
70 per cent of the halo effective turn-around time. This makes it a useful qualitative 
tool for understanding certain average properties of galaxies, such as total spin and 
angular-momentum distribution within haloes, but with a random scatter of the order 
of the signal itself. Non-linear changes in spin direction cause a mean error of ~ 50° 
in the TTT prediction at t^, such that the linear spatial correlations of spins on 
scales > 1 /i~^Mpc are significantly weakened by non-linear effects. This questions the 
usefulness of TTT for predicting intrinsic alignments in the context of gravitational 
lensing. We find that the standard approximations made in TTT, including a second 
order expansion of the Zel'dovich potential and a smoothing of the tidal field, provide 
close-to-optimal results. 

Key Vifords: cosmology: dark-matter - cosmology: large-scale structure of Universe 
- cosmology: theory - galaxies: formation - galaxies: haloes - galaxies: structure 



1 INTRODUCTION 

Angular momentum is clearly one of the key physical ingre- 
' dients in the process of galaxy formation. The total angular 
momentum, and its distribution, must have a crucial role in 
determining the galaxy history and final type. It has there- 
fore been a subject for classical investigations, pioneered by 
Hoyle (1949), and then analysed qualitatively by Peebles 
(1969). This led to the 'standard' theory for the origin of 
angular momentum in the framework of hierarchical cosmo- 
logical structure formation, the tidal-torque theory (TTT) 
due to Doroshkevich (1970) and White (1984). 

Special interest in the subject has been revived recently 
because of a "spin crisis" , arising from cosmological simula- 
tions of galaxy formation which use hydrodynamical gravita- 
tional codes to follow the gas dynamics inside dark-matter 
haloes, and semi-analytical recipes for star formation and 
feedback. These simulations seem to yield luminous galax- 
ies that are significantly smaller, and of much less angular 
momentum than observed disc galaxies (Navarro, Frenk & 
White 1995; Navarro & Steinmetz 1997, 2000). The discov- 
ery of massive black holes in galactic centers (Kormendy 
& Richstone 1995; Magorrian et al. 1998; Gebhardt et al. 



2000a,b) provides another motivation for understanding the 
detailed distribution of angular momentum within galax- 
ies. A general need for a detailed recipe for the build-up of 
angular momentum in galaxies comes from the developing 
semi-analytic models, which have been proven very useful 
in the study of galaxy formation (e.g. White & Frenk 1991; 
Kauffmann, White & Guiderdoni 1993; Cole et al. 1994; 
Somerville & Primack 1999). Another current motivation 
comes from weak lensing studies. Our ability to reconstruct 
maps of the cosmic mass distribution from weak gravita- 
tional lensing measurements (Bacon, Refregier & Ellis 2000; 
Kaiser, Wilson & Luppino 2000; van Waerbeke et al. 2000, 
2001; Wittman et al. 2000) is hampered by the unknown 
intrinsic distribution and alignment of galaxy shapes. As a 
first guess, one crudely assumes that the shapes and ori- 
entations of background galaxies are uncorrelated in space, 
and that the detected correlations of galaxy ellipticities are 
solely due to the lensing by the foreground mass distribu- 
tion. However, a series of observations (Brown et al. 2000; 
Pen, Lee & Seljak 2000) and numerical simulations (Croft 
& Metzler 2000; Heavens, Refregier & Heymans 2000; Dekel 
et al. 2000) evidenced that both the intrinsic galaxy shapes 
and their spins are spatially correlated. Theoretical models 
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for the alignment of disc galaxies orientations tend to as- 
sume that they are induced by spin correlations (Catelan, 
Kamionkowski & Blandford 2001; Crittenden et al. 2001). 

These issues add a timely aspect to the motivation for 
revisiting the classical problem of angular momentum, in an 
attempt to sharpen our understanding of its various com- 
ponents. This includes the first step of angular-momentum 
acquisition by dark haloes, which might be particularly as- 
sumed to be fairly well understood. The basic idea of TTT 
is that most of the angular momentum is being gained grad- 
ually by proto-haloes in the linear regime of density fiuctu- 
ations growth, due to tidal torques from neighboring fiuctu- 
ations, and that this process continues until the proto-halo 
reaches its maximum extent. The further assumption is that 
only little angular momentum is being exchanged between 
haloes later on in the non-linear regime, after the proto- 
haloes have decoupled from the expanding background and 
collapsed to virialized systems. It is commonly assumed that 
the baryonic material, which in general follows the dark- 
matter distribution inside each proto-halo, gains a similar 
specific angular momentum and carries it along when it con- 
tracts to form a luminous galaxy at the halo centre. This 
should allow us to predict galactic spins using the approx- 
imate but powerful analytic tools of quasi-linear theory of 
gravitational instability. 

In a series of papers, we evaluate the performance of the 
TTT approximation, and trace the roles of its various ingre- 
dients, using a cosmological A^-body simulation with thou- 
sands of well-resolved haloes. As we do so, we find to our 
surprise that some of the basic ingredients of TTT, which 
are commonly assumed to be of 'text-book' status, involve 
certain unjustified assumptions and confused understand- 
ings, which may lead to poor approximations. These papers 
represent attempts to clarify some of these controversial is- 
sues. 

We address the TTT at different levels. In this paper 
(Paper I), we evaluate how well does the approximation pre- 
dict the final angular momentum of a halo, given full knowl- 
edge of the corresponding initial proto-halo and the cosmo- 
logical realization. In Paper II (Porciani, Dekel & Hofi'man 
2001), we attempt a deeper level of understanding of the 
origin of halo angular momentum, by investigating the rela- 
tion between the different components of TTT. This study 
connects to the fundamental open question of how to iden- 
tify a proto-halo in a given realization of initial conditions. 
In Paper III (Porciani & Dekel in preparation), we revise 
the standard scaling relation of TTT, based on the fact that 
the shear tensor is only weakly correlated with the density 
contrast of the proto-halo (the density and shear being the 
trace and traceless parts of the same, deformation tensor). 
This scaling relation is used to predict the typical angular- 
momentum profile of haloes (Dekel et al. 2001; Bullock et 
al. 2001b), and to provide a simple way to incorporate angu- 
lar momentum in semi-analytic models of galaxy formation 
(Mailer, Dekel & Somerville 2002). 

Given the initial conditions of a proto-halo and its en- 
vironment, TTT predicts the final halo angular momentum 
based on four specific assumptions as follows: (a) The fiow 
is laminar, with a one-to-one correspondence of Eulerian 
and Lagrangian positions, or the velocity field is properly 
smoothed such that laminarity prevails in practice, (b) The 
velocities obey the Zel'dovich approximation, (c) The po- 



tential at every point within the proto-halo can be approxi- 
mated by its Taylor expansion to second order (with respect 
to spatial separation) about the centre of mass, (d) There 
is little contribution to the halo angular momentum from 
non-linear effects. In the current paper, we set to evaluate 
the global success of TTT in this task, and try to address 
the validity of each of the above assumptions. The role of 
non-linear efi'ects, in particular, turns out to be the key is- 
sue. This will be an attempt to follow and improve on earlier 
work by White (1984), Hoffman (1986, 1988), Barnes & Efs- 
tathiou (1987), Catelan & Theuns (1996), Lee & Pen (2000) 
and Sugerman, Summers & Kamionkowski (2000). 

The outline of this paper is as follows: In §d we summa- 
rize the basics of linear tidal-torque theory. In § J we describe 
the simulation and the halo finder. In §^ we describe the im- 
plementation of TTT to proto-haloes. In §P we evaluate the 
success of TTT in predicting the amplitudes of halo spins. In 
^ we assess the success of TTT in predicting the directions 
of halo spins. In ^ we address the spatial coherence of spin 
directions, which is directly relevant to intrinsic alignments 
in the context of gravitational lensing. In §^ we discuss our 
results and conclude. 



2 TIDAL-TORQUE THEORY 

We start by summarizing the basics of linear tidal- 
torque theory. The framework is the standard Friedmann- 
Robertson- Walker cosmological model in the matter- 
dominated post-recombination era, in which the matter is 
assumed to be a collision-less fluid and the dynamics is 
driven by gravity only. In the linear regime, the developing 
proto-haloes are assumed to be small perturbations about 
the mean universal density. The peculiar velocity field, due 
to the growing mode of gravitational instability, describes a 
potential flow with no curl, assuming that any primordial 
source of vorticity (e.g., due to some primordial turbulence, 
see Jones 1976) has decayed away due to the expansion of 
the universe. 

Given a proto-halo, a patch of matter occupying an Eu- 
lerian volume 7 that is destined to end up in a virialized 
halo, the goal is to compute the halo angular momentum 
about the centre of mass, to the lowest non- vanishing order 
in perturbation theory. The angular momentum at time t is 
deflned by 



L(t)= / p{r,t)[r{t) 



,(t)]dV, (1) 



where r and v are the position and peculiar velocity vectors, 
and the centre-of-mass quantities, rem and Vcm, are defined 
as usual. The term proportional to Vcm does not contribute 
to L, and it will not be considered any further. Re-write 
Eq. (^) in comoving units, x = r/a(t) and v — adx./dt, 
where a{t) is the universal expansion factor, and replace 
p{x,t) with the density contrast 5(x,t) = p{-x.,t) / pit) — 1 
with respect to the average density pit), to obtain 

L(t) =p(t)a^(t) / [l-h5(x,t)] [x(t) -Xcm(t)] X xd^x. 

J -y 

(2) 

A dot denotes a derivative with respect to cosmic time, t. In 
the matter-dominated era, p{t)a?it) = poaf, = const, where 
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the subscript denotes quantities evaluated at the present 
time. 

In the Lagrangian description of fluid dynamics, the 
comoving Eulerian position of each fluid element is given 
by its initial, Lagrangian position q plus a displacement: 
x(q, t) = q + S(q, t). When fluctuations are small, or when 
the flow is properly smoothed, the mapping q ^ x is re- 
versible and the flow is laminar. Then the Jacobian deter- 
minant J — ||9x/9q|| does not vanish, and the continuity 
equation implies 1 -I- 5[x(q, t)] = J^^(q, t). Substituting in 
Eq. (0) we obtain, for a laminar flow, 

L(t) = a\t) poal [ [q-q + S(q,t)-Sl xS(q,t)d^g, 



and the inertia tensor of F, 



(3) 

where F is the region in Lagrangian space corresponding to 
7 in Eulerian space. Barred quantities are averages over q in 
F. Note that Eq. is exact in the absence of orbit crossing. 

The displacement S is now spelled out using the 
ZeVdovich approximation (Zel'dovich 1970), which is lin- 
ear in Lagrangian terms and mildly non-linear in Eulerian 
space. It assumes that the proportionality between the ve- 
locity potential and the gravitational potential 0(q, t), which 
holds in the linear regime for the growing mode of density 
fluctuations 5 oc D{t), can be extended into the mildly non- 
linear regime. This implies that the spatial and temporal 
dependences of S can be decoupled, 

S(q,t) = -D(f)V<l>(q). (4) 

where $(q) = 0(q, t)/47rGp(t)a^(t)D(t) (with G Newton's 
gravitational constant). It also implies that the peculiar 
velocity and acceleration fields are parallel. Substituting 
Eq. (M) into Eq. (H) one obtains 



L(f) = ~a\t)D(t) poal / (q - q) x V-l'(q) d?q 



(5) 



We see that the explicit growth rate is L oc a^{t)D{t). For 
an Einstein-de Sitter universe this is oc f.Q 

Next, assume that the potential is smoothly varying 
within the volume F, such that it can be approximated by 
its second order Taylor expansion about the centre of mass 

q, 



<f.(q') ~ $(0) + — 



1 



2 dq'.dq'^ 



(6) 



q'=0 



where q' = q — qj^ Substituting in Eq. (|^) one obtains the 
basic TTT expression for the ith Cartesian component: 



Li{t) = a {t)b{t) Eijk T>ji Ilk , 



(7) 



where eijk is the fully antisymmetric rank-three tensor, and 
the two key quantities are the deformation tensor at q' = 0, 



dq'idq'. 



(8) 



polo / q'i q'j d^q 



(9) 



Note that only the traceless parts of the two tensors mat- 
ter for the cross product in Eq. (17|). These are the veloc- 



and the 
(/«/3)5y. Thus, 



ity shear or tidal tensor, Tij = Vi 
traceless quadrupolar inertia tensor, li- 
to the first non- vanishing order in perturbation theory, angu- 
lar momentum is transferred to the developing proto-halo by 
the gravitational coupling of the quadrupole moment of its 
mass distribution with the tidal field exerted by neighboring 
density fluctuations. The torque depends on the proto-halo 
shape, the external tidal field, and the misalignment be- 
tween the two. In particular, L vanishes if F is a sphere or 
is bounded by an equipotential surface of $. 

There is no a priori justification for truncating the ex- 
pansion of the potential after second-order, Eq. (n). More 
generally, the expansion includes further terms. 



U{t) 



with 



(3)/ 



(10) 



a\t)D{t)e,,k 



ir)(3) r(3) , 
2 jltn link ~^ ' 



dq'idq'.dq'k 



r(3) - 3 



q'=0 



q^q 



(11) 



Here, multipole moments of a given order of the proto-halo 
mass distribution are coupled with spatial derivatives of the 
same order of the potential. This is a consequence of the 
scaling of gravitational interactions, and indeed, the 
total torque for a localized charge in an external electrostatic 
field admits a similar expansion (e.g., Jackson 1975, exercise 
4.5). Note that dipole interactions [i.e., terms proportional 
to p X g(q' = 0) with p the mass dipole moment and g the 
gravitational acceleration] do not contribute to the torque 
because the spin is considered in the centre-of-mass frame, 
where the mass dipole vanishes. On the other hand, dipole 
interactions are expected to play an important role in the 
evolution of orbital angular momentum. 

After the dark-matter haloes have detached from the ex- 
panding background and collapsed into condensed systems, 
the tidal torquing continues, but at a highly reduced rate. 
Peebles (1969) estimated the typical torque exerted on a 
quadrupolar mass distribution by a population of point-like 
galaxies, showing that the transfer of angular momentum be- 
comes less and less efficient with time. This is a consequence 
of the fact that haloes become small when they collapse, as 
well as the fact that they tend to move away from each other 
due to the Hubble expansion. 



q'=0 



* A numerical integration of the cosmological relations shows 
that for all flat cosmologies the approximation a^D oc D^^^ holds 
to an excellent accuracy, and it deviates only slightly for open 
cosmologies. 

t This is equivalent to assuming that the velocity field is well 
described by its linear Taylor expansion, i.e., Vi — Vi ~ f ij 9j 
(where v is in comoving units). 



3 SIMULATION AND HALO FINDING 

The A''-body simulation analysed here was performed as part 
of the GIF project (e.g., Kauffmann et al. 1999) using the 
parallel, adaptive particle-particle/particle-mesh (AP'^M) 
code developed by the Virgo consortium (Pearce et al. 1995; 
Pearce & Couchman 1997). As an example, we use a simu- 
lation of the rCDM scenario, in which the cosmology is flat. 
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Einstein-de Sitter with density parameter fim = 1, and Hub- 
ble constant Ho = 100 /i kms"^ Mpc"^ with h = 0.5. The 
power spectrum of initial density fluctuations is CDM with 
shape parameter F = 0.21. The amplitude of the power spec- 
trum is set such that the rms overdensity in top-hat spheres 
of radius 8/i~^Mpc, as extrapolated using linear theory to 

the present epoch, is erg = 0.51.^ 

The simulation was performed in a periodic cubic box of 
side 84.55 /i~^Mpc, with the mass represented by 256"^ par- 
ticles of 1.0 X 10^°h~'^MQ each. Lon B-ranee eravitational 
forces were computed on a 512^ mesh, while short-range 
interactions were calculated as in Efstathiou & Eastwood 
(1981). At late times {z ^ 3), the corresponding gravita- 
tional potential asymptotically matches a Plummer law with 
softening e = 36h~^kpc both at large and small scales. The 
simulation started at a redshift z — 50 and ended at the 
present time, z — 0. The initial conditions were generated 
by displacing particles according to the Zel'dovich approxi- 
mation from an initial stable 'glass' state (e.g. White 1996). 
More details are given in Jenkins et al. (1998). 

We identify virialized dark-matter haloes at 2: = as 
follows. First, we select an initial set of groups of particles 
using a friends-of-friends method (e.g., Davis et al. 1985) 
with a linking length b = 0.2 in units of the average inter- 
particle distance. This algorithm identifies regions bounded 
by a surface of approximately constant density contrast, 
p/p ~ 3/(27rfe'') (e.g. Cole & Lacey 1996). Assuming that the 
density profiles of these groups approximate singular isother- 
mal spheres, the selection criterion can be expressed in terms 
of the average density contrast in the halo, p/p ~ 9/(27rb'^). 
Thus, b = 0.2 corresponds to haloes with a mean density 
contrast ~ 179, in general agreement with the predictions of 
the spherical collapse model for perturbations whose outer 
shells have collapsed recently. 

We then remove unbound particles from each halo by 
the following iterative procedure. At each step, we compute 
the energies of all particles in the current centre-of-mass 
frame, and remove the most unbound particle. We then re- 
compute the energies of all the particles and so on. The iter- 
ation stops when no unbound particles are found. Note that 
a particle that has been removed is allowed to re-enter the 
halo if it becomes bound at a later step due to changes in the 
position and velocity of the centre of mass. Clumps which 
end up containing between 100 and 200 particles loose, on 
average, ~ 5.5 per cent of their mass as unbound parti- 
cles. In a few cases, typical of small haloes, the friends-of- 
friends algorithm identifies 'fake' haloes which are not bound 
gravitationally. For example, seven friends-of-friends haloes 
(out of ~ 7300) composed by more than 100 particles end 
up containing only 10 bound particles or less after the un- 
binding procedure. More massive haloes are, naturally, less 

t Such a power spectrum may arise, for example, if a late de- 
cay of massive r neutrinos produces an additional background of 
relativistic e and fi neutrinos, thus delaying the onset of matter 
domination and lowering the effective value of T (Bond & Efs- 
tathiou 1991; Efstathiou, Bond & White 1992; White, Gelmini 
& Silk 1995). Here it simply serves us as a scenario which obeys 
many of the observational constraints from large-scale structure, 
including the normalizations set by cluster abundance and by 
COBE's measurements of large-angle anisotropics in the CMB, 
independently of its physical origin. 



affected by unbound interlopers because their potential wells 
are typically deeper. For instance, the mean fraction of un- 
bound mass removed from clumps which end up containing 
between 200 and 1000 (1000-3000) particles is ~ 2.9 per cent 
(~ 2.0 per cent). In the following we consider only haloes 
which contain more than 100 particles after the un-binding 
procedure is applied. Our conclusions regarding TTT turn 
out to be insensitive to the removal of unbound particles, 
and to increasing the minimum-mass threshold. In order to 
try to restrict ourselves to galactic haloes, we also exclude 
all the ~ 200 haloes which consist of more than 3000 par- 
ticles - these tend to have more substructure and are likely 
to correspond to groups or clusters of galaxies. 

The robustness of our conclusions regarding TTT with 
respect to the halo-finding algorithm (e.g., friends-of-friends 
versus fitting a spherical or ellipsoidal density profile. Bul- 
lock et al. 2001a) is investigated in another paper in prepa- 
ration. Note that all the haloes in our current sample are not 
subclumps of larger host haloes, for which linear theory is 
not expected to be valid. This excludes about 10 per cent of 
the haloes that are more massive than lO^'^ h~^M^ (Sigad 
et aL 2001). 



4 IMPLEMENTING TTT 

The proto-halo regions F in Lagrangian space are defined 
straightforwardly by tracing all the halo particles into their 
Lagrangian positions. Generally, most of the proto-halo mass 
is contained in a simply connected Lagrangian region, but in 
some cases the proto-halo may be divided into several com- 
pact regions which are connected by thin filaments. Nearly 
10 per cent of the proto-haloes are characterized by extended 
filaments departing from a compact core. 

We stress that the proto-haloes at earlier times are de- 
fined solely based on their identification as virialized haloes 
at z = 0, without worrying about the details of the merger 
history, and without applying any other criteria at the initial 
conditions. Some of the results reported below are specific 
to this choice of considering all today's halo particles at all 
times rather than, e.g., investigating the history of the halo 
major progenitor (Vivitska et al. 2001; Wechsler et al., in 
preparation). 

For each proto-halo, we compute the Lagrangian inertia 
tensor by direct summation over its A'' particles of mass m 
each: 

h,=mf2l[^"^l}"\ (12) 

n=l 

with q the position of the n-th particle with respect to 
the halo centre of mass. 

We use three alternative methods to measure the shear 
tensor at the proto-halo centre of mass in the initial condi- 
tions of the A-body simulation. 

Method 1, top-hat smoothing: We smooth the gravita- 
tional potential used to generate the initial Zel'dovich dis- 
placements, and then differentiate it two times with respect 
to the spatial coordinates to obtain the deformation tensor. 
Smoothing is done using a top-hat window function, while 
derivatives are computed on a grid. The top-hat smooth- 
ing radius in comoving units, 7?th, is taken to be defined by 
(47r/3)poi?th ~ M, with M the virial halo mass. This is the 
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only practical method for application in analytic and semi- 
analytic modelling of galaxy formation. We refer to is as the 
standard method of TTT. 

Method 2, fit total: The smoothing applied in Method 
1 is associated with a loss of information. First, the coarse- 
graining procedure averages over small-scale structure. Sec- 
ond, the spherical smoothing kernel does not maintain the 
details of the proto-halo shape, since, as we find in Paper II, 
most proto-haloes are far from being spherically symmetric. 
We minimize smoothing in the computation of the tidal ten- 
sor in the following way. The comoving linear velocity field, 
in the proto-halo centre-of-mass frame, is expanded in terms 
of the corresponding Eulerian comoving separation, x, as 

Vi — T>ij Xj . (13) 

The components of the deformation tensor Dij = dvi/dxj 
are computed by least-squares fitting to the simulation data 
(cf. Kaiser 1991). This is realized by solving the system of 
equations 



(")^(n) 

i 1 



Ti — 1 n— 1 

5][©,,4"'4")-fo,,4"'x("'] = (14) 

n = l 

JV 
n = l 

Since the linear growing mode of the gravitational velocity 
field is curl-free, Eq. (114) has been derived after imposing 
the condition Dij — Dji. We then compute the symmetric, 
traceless shear tensor Tij = Dij — {l/3)DiiSij. The shear 
field computed this way corresponds to minimal smooth- 
ing, at the level of the individual particles in the A'^-body 
simulation. However, even though no velocity averaging is 
performed, method 2 is, in some sense, similar to using 
a smoothing kernel which adapts itself to the proto-halo 
shape. The resulting Tij describes the quadrupolar structure 
of the velocity field in the proto-halo at the particle level. 
Method 2 is therefore expected to provide better predictions 
for L than Method 1 whenever the top-hat smoothing is as- 
sociated with a severe loss of information. 

Method 3, fit external: The interactions between par- 
ticles that belong to the same proto-halo cannot create 
any net angular momentum. This motivates us to consider 
the shear tensor generated solely by the density perturba- 
tions lying outside the proto-halo volume, F. This is real- 
ized by decomposing the linear velocity field into two parts. 



/self) 



foxt) 



accounting for the gravitational pull of 
the mass lying in and outside the proto-halo, respectively. 
In practice, for each particle in F, we first compute the lin- 
ear gravitational acceleration generated at 2 = 50 by all the 
other point masses forming the proto-halo. ^From this quan- 
tity, we then subtract the acceleration caused by the same 
particles as if they were part of a uniform background. This 
is computed using the unperturbed particle positions of the 
glass state which, formally, correspond to 2: — > 00. Rescaling 
the result we obtain v'^"'^'. The difference between the ac- 
tual velocity of the particles in the simulation at 2 = 50 and 
this "self gravitationally induced term is taken to represent 
.^(oxt) jjj(jgg{j^ estimate for v'™*-* generates, in practice, 
the same spin as the total velocity field. Eventually, the ex- 
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Figure 1. Predicted versus actual angular-momentum amplitude 
at the initial time, z = 50. Shown is the distribution over all 
the haloes in our sample. The three histograms correspond to 
the different methods of computing the shear tensor discussed 
in the main text (§^: using top-hat smoothing (thick line) or 
fitting with no smoothing of either the total shear (thin line) or 
the external part only (dotted line). To improve readability, the 
thin and dotted histograms have been slightly displaced in the 
horizontal direction. 



ternal tidal field is computed by least-squares fitting. Note 
that, since 2?'^""'' oc lij, the self field does not contribute to 
the angular momentum also at a linear level - i.e. there is 
no self-torquing. 

The results obtained by Method 2 and Method 3 con- 
tain spurious contributions due to errors in the determina- 
tion of the deformation tensor by the fit of a linear model to 
a discrete set of points. However, since the external shear is 
expected to vary more smoothly than the self gravity on a 
proto-halo scale, should be described more accurately 

by the linear expansion of the velocity field. For this reason. 
Method 3 should provide a better estimate of the angular 
momentum. Note, however, that additional noise is intro- 
duced in Method 3 by replacing the uniform background 
density distribution with the glass state. 



5 SPIN AMPLITUDE 

Given a proto-halo at the initial conditions, we first investi- 
gate the performance of TTT, Eq. (^), in terms of predicting 
the amplitude of the halo angular momentum. Any system- 
atic deviation of the mean predicted amplitude from that 
obtained in the simulations can be eliminated by adjusting 
the effective time near turn-around at which the linear TTT 
growth is assumed to stop. Random deviations, both due to 
the approximations made already at the initial conditions as 
well as non-linear effects at late epochs, cannot be avoided 
and should be estimated. 

We start by testing the accuracy of TTT in the lin- 
ear regime, at the initial time, z = 50. For each halo, we 
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compute the inertia tensor lij and the shear tensor Tij us- 
ing the three different methods described in §^. We then 
compute the amplitude Lmodci as predicted by Eq. (^), and 
compare it with the actual angular momentum of the proto- 
halo particles, Ltruc- The distribution of the ratio of the two 
over the protohaloes is shown in Fig. ^, and characterized 
in Table |l|. The mean (median) values for methods 1, 2, 3 
are 0.95 (0.83), 1.05 (0.99), 1.01 (0.95) respectively. In all 
cases the probability distribution is positively skewed, and 
is leptokurtic (with a sharp peak and long tails) with re- 
spect to a Gaussian distribution in log-space. It is therefore 
not very useful to quantify the dispersion about the mean 
by a single standard deviation, so we quote instead the 68.3 
per cent confidence interval, obtained by locating the values 
corresponding to the 15.9 and 84.2 percentiles in the cumu- 
lative distribution. For methods 1, 2, 3 we find (0.53, 1.28), 
(0.72, 1.31), (0.71, 1.24) respectively. As a single measure of 
scatter we use half the 68.3 per cent range, which is 0.37, 
0.29, 0.26 respectively. All three methods involve the second- 
order truncation of the expansion of the potential and some 
smoothing. These are the only sources of error, since the 
Zel'dovich approximation is exact at the initial conditions 
of the simulation by construction. Even though the cases 
with minimum smoothing (methods 2 and 3) involve noisy 
fits of the small-scale structure, they give somewhat more 
accurate predictions. The external fit does just a little bet- 
ter than the total fit, both in terms of systematic error and 
scatter. However, the smoothed tidal field (method 1), which 
is the only practical method for semi-analytic modelling of 
galaxy formation, also leads to a small systematic deviation 
of the mean, of only ~ 5 per cent. The scatter, though, is of 
the order of the signal for all three methods already at the 
initial conditions. 

These numbers change only slightly when restricting 
the analysis to the very well-resolved haloes containing 1000 
particles or more. In this case, the mean (median) of the dis- 
tribution of the amplitude ratio is 0.90 (0.79), 1.06 (1.02), 
1.02 (0.98) for methods 1, 2, 3 respectively. The correspond- 
ing scatter is 0.32, 0.30, 0.27. 

When we add the third-order term in the Taylor ex- 
pansion of the gravitational potential, i.e., the second term 
in the right-hand side of Eq. (^l|), we find no significant 
improvement in the systematic error and in the scatter. In 
this case the mean (median) and scatter of Lmodei/itruc are 
0.95 (0.83) ±0.36, almost identical to the second-order re- 
sults (top-hat smoothing is adopted in both cases). A more 
detailed inspection reveals the following. When consider- 
ing only those haloes for which standard second-order TTT 
over-predicts or under-predicts the spin amplitude by more 
than a factor of 3, we find that by adding the third-order 
term one improves the prediction in about 75 per cent of the 
cases, though the overall scatter increases. In the case of se- 
vere over-prediction, the mean (median) and scatter are 4.54 
(3.41) ±0.44 for second-order TTT and 3.81 (2.63) ±0.79 
for third-order TTT. In the case of severe under-prediction, 
the numbers are 0.25 (0.26) ±0.07 for second-order TTT 
and 0.40 (0.33) ±0.16 for third-order TTT. On the other 
hand, when the second-order TTT predictions are reason- 
ably good, the addition of the third-order term makes no 
general improvement. For the haloes were the prediction is 
within 30 per cent of the true value, the mean (median) and 
scatter are 0.94 (0.92) ±0.17 for second-order TTT and 0.94 
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Figure 2. Probability distributions for the angle (top) and the ra- 
tio of the amplitudes (bottom) between the halo spin predicted by 
TTT considering the first term on the right hand side of Eq. ( |ll| ) , 
L(^', and the second term, L^'^). The dashed line in the top panel 
is the expected distribution for random orientations. 



(0.91) ±0.21 for third-order TTT. Part of the failure of the 
third order term to improve the results in these cases could 
be attributed to numerical errors in the computation of the 
third-order term, which are particularly important when the 
high-order corrections are small. 

Fig. shows the distribution of the angle between the 
direction of L as predicted by standard second-order TTT 
(top-hat smoothing) and by the third-order term alone. We 
see that the additional contribution points at an almost ran- 
dom direction compared to the leading-order term, with an 
average angle of 98° and a scatter of ~ 45° . Fig. |^ shows the 
distribution of the corresponding amplitude ratio, L'^'/L'^', 
with an average of 0.39 and a scatter of 0.27. These indicate 
that any truncation of the expansion is doomed to intro- 
duce an error that is not negligible compared to the signal. 
We conclude that the second-order expansion is the most 
sensible and practical procedure, but it is associated with a 
noticable error. 

In the second step of our testing, we address the growth 
rate indicated by TTT in the quasi-linear regime, L oc a^D. 
We do it by comparing for each proto-halo the model and 
true L values at z = 3, where the fluctuations are still linear 
on the Lagrangian scale of the protohaloes, but the proto- 
halo substructures, and maybe even the haloes themselves, 
already involve non-linear fluctuations. We show in Fig. ^ 
the distribution of the ratio Lmodoi to Ltrue for standard 
TTT. We show, for comparison, the corresponding distri- 
bution for the case where the model angular momentum is 
computed directly from the z = 3 velocities, as predicted 
by the Zel'dovich approximation from the initial conditions 
with the corresponding growth. 

The deviation of the ZA histogram from a Dirac S func- 
tion at unity reflects the deviation of the angular-momentum 
growth rate from the linear prediction, L oc a^D. It does 
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Figure 3. Model versus true angular-momentum amplitudes in 
the quasi-linear regime, z = 3. Shown is the distribution over 
all the haloes in our sample. The thick histogram corresponds to 
standard TTT based on the truncated and smoothed shear tensor. 
The thin histogram (which has been slightly shifted to the right) 
corresponds to L computed directly from the z = 3 Zel'dovich 
velocities. The deviations from unity in the latter is solely due to 
the non-linear effects in the growth rate of angular momentum. 



not involve the errors associated with the smoothing and 
the truncation of the potential because the initial condi- 
tions also obey the Zel'dovich approximation. The predicted 
growth rate is slightly higher than the true growth rate, but 
the corresponding 8 per cent systematic overestimate and 
the 0.19 scatter are relatively small. The additional effects 
of truncation plus smoothing in TTT are responsible for 
the larger scatter in the TTT histogram. In this case, the 
systematic error is only 2 per cent, due to a slight system- 
atic underestimate by TTT of the spin amplitude at z = 50, 
which partly compensates for the overestimate of the growth 
rate. Note, however, that the median of the distribution (see 
Table |l|) keeps roughly constant between z = 50 and z — 3. 
The scatter of 0.43 is comparable to the signal, and is slightly 
larger than the initial scatter of 0.38 due to the deviations 
from the linear growth rate. 

In the third step, we compare the spin amplitudes at z = 
0, after non-linear processes have affected the true spin. We 
first wish to determine the effective time in each halo history, 
or on average for all haloes, at which the linear growth of 
angular momentum a la TTT should be assumed to stop 
in order to obtain the best prediction for its amplitude at 
z = 0. Then, we wish to evaluate the scatter about the mean 
TTT resuh. 

In order to learn about the actual evolution in the non- 
linear regime, we show in Fig. ^ the distribution over the 
haloes of the ratio between Ltruo at 2 = (t = to) and at 
z = 3 {t = to/S). The actual growth factor between these 
times has a mean value of 2.9 (with a scatter of ^ 1.4), which 
implies that the linear growth should be stopped at an effec- 
tive time 2.9 to/8 — 0.36 to (or z ~ 1). This time is somewhat 



Figure 4. True growth factor of angular momentum between 
^ = 3 and 2 = 0. Shown is the distribution over all the haloes in 
our sample. 

earlier than maximum expansion for the outer shell of the 
halo (using the terminology of the spherical collapse model) , 
because our haloes, by definition, have their outermost shells 
reaching virialization today, and were therefore at maximum 
expansion at t ~ to/2 (or z ~ 0.6). Linear growth until max- 
imum expansion would have implied a growth factor of 4 
instead of 2.9 since z — 3. 

The actual evolution of spin amplitudes of the simu- 
lated haloes is shown in Fig. |5[ The average and scatter 
refer alternatively to the ratio of L(t) and the initial proto- 
halo spin, and the ratio of L{t) and the final halo spin. We 
see that, in agreement with the results shown in Fig. ^ the 
best effective time for stopping the TTT growth is ~ 0.3 to. 
This would properly adjust the mean, but would leave a sig- 
nificant scatter from halo to halo as indicated by the error 
bars. 

In Fig. 1^, we show the distribution of the ratio of model 
to true L values at z = 0, with the TTT growth stopped in 
three different ways. In the top panel, the TTT growth of 
each halo has been stopped at a fixed fraction of its turn- 
around time tta, as predicted by the spherical collapse model 
based on the value of the initial density fluctuation smoothed 
on the scale of the proto-halo (namely, the time at which the 
linearly extrapolated density contrast is 5 = 1.06). By stop- 
ping the TTT growth for each halo at tTTT ~ 0.56 tta, the 
median of the ratio l/modoi/itruc becomes unity. Alterna- 
tively, in the central panel, the TTT growth for each halo 
has been stopped at a fixed fraction of the time tiD when 
the smoothed proto-halo is predicted by the Zel'dovich ap- 
proximation at the initial conditions to collapse along the 
first principal axis of the deformation tensor. Requiring that 
the median of the corresponding amplitude ratio is equal to 
unity gives tTTT = 0.39 tiD. In the bottom panel, the lin- 
ear spin growth has been stopped at the same epoch for 
all the haloes, tTTT = 0.35 to corresponding to z = 1.03, 
which has been determined again to bring the median value 
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Table 1. Characteristics of the distribution of the ratio of amphtudes of the spins as predicted 
by the model and as measured in the simulations. Listed are the mean and several percentiles, 
as well as the scatter as measured by half the central 68.3 per cent region. TTT=tidal-torque 
theory, NB=A'^-body simulations, ZA=Zerdovich approximation. 

Comparison z Fig. Mean 2.5% 15.9% 50% 84.2% 97.5% Scatter 
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Figure 5. True evolution of spin amplitude in the simulation 
compared to TTT. Shown are the average and the 68.3 per cent 
confidence interval over the haloes at different times, of the ratio 
of spin amplitude L{t) and either the spin today (top panel) or the 
spin at the initial conditions (bottom panel) . The continuous lines 
represent the linear TTT growth stopped at maximum expansion 
{t = 0.5 to). The dotted lines mark the average value at tg. 



of imodci/^truc to uiilty. As at z = 3, we show the his- 
tograms for TTT and for the pure Zel'dovich velocities at 
tTTT- The corresponding percentiles are listed in Table |l[ 
Note that, despite the different stopping times, the three 
probability distributions are very similar; the 68.3 per cent 



confidence level corresponds to a factor of 2 random error. 
The histograms for TTT and for the pure Zel'dovich veloci- 
ties show a similar scatter, indicating that the choice of Ittt 
is the dominant factor in determining the scatter compared 
to the other approximations involved in TTT. 

The impression from Fig. ^ that the growth of spin 
amplitude really stops at some specific epoch is somewhat 
misleading. When we investigate halo by halo, we find that 
about 38 per cent of our haloes actually loose angular mo- 
mentum between z = 0.6, z = 0.2 and z = due to non- 
linear interactions (see a similar finding by Sugerman et al. 
2000). On average, these haloes loose ~ 30 per cent of their 
spin between z = 0.6 and z — 0. This is partly a result 
of considering in each proto-halo the same group of parti- 
cles at all times, based on their association with a halo at 
z = 0. Angular momentum from these proto-halo particles 
may be transferred to outer shells of matter. Mergers, tidal 
encounters and the presence of substructure all contribute to 
angular-momentum transfer. This decrease of L with time 
at late epochs indicates that in many case the TTT growth 
actually continued roughly until turn-around, stopped, and 
started declining thereafter. This is one of the reasons for 
why the efi'ective time for stopping the TTT growth should 
be somewhat before turn-around. 

Some fraction of the haloes have their spin growing also 
at late epochs: about 21 per cent of our spins grow between 
z — 0.6, z — 0.2 and z — 0. On average, these haloes in- 
crease their spin by ~ 60 per cent between z = 0.6 and 
z — 0. However, for about 57 per cent of our haloes the evo- 
lution of spin amplitude in the non-linear regime is weak, 
less than 20 per cent between z — 0.6 and z = 0, in pleas- 
ant agreement with the concept of TTT. Considering the 
whole halo population, we find that the mean value and the 
scatter of the ratio L{z = 0.6)/L{z = 0) are 1.19 and 0.39 
- or 1.01 and 0.34 for the inverse ratio. We conclude that 



© 0000 RAS, MNRAS 000, 000-000 




Figure 6. Model versus true spin amplitudes in the non-linear 
regime, at z = 0. The TTT growth has been stopped at txTT 
such that the median of the ratio imodel/^true is unity. In the 
top panel fxTT =0.56 tta, a fixed fraction of the turn-around time 
as predicted for each halo by the spherical-collapse model based 
on the smoothed initial density fluctuation. In the middle panel 
^TTT = 0.39 tiD, a fixed fraction of the time of orbit crossing 
along the first principal axis of the deformation tensor as pre- 
dicted by the Zel'dovich approximation based on the smoothed 
initial conditions. In the bottom panel, the linear spin growth has 
been stopped at the same epoch for all haloes, txTT = 0.35 to. 
The thick histogram corresponds to standard TTT, with a scatter 
of a factor of two. The thin histogram (which has been slightly 
displaced to the right) corresponds to Lmodol computed directly 
from the Zel'dovich velocities at ixTT! it does not involve the 
approximations made in TTT at the initial conditions and can 
be interpreted as reflecting only the scatter due to the choice of 
txTT • 



even though non-linear effects do not significantly change 
the average amplitude of halo spins, the variety of evolu- 
tionary paths result in a large scatter about the mean TTT 
prediction at 2 = 0. 

The bottom line is that standard TTT provides a suc- 
cessful order-of-magnitude estimate for the halo angular- 
momentum amplitude. If TTT growth is stopped at turn- 
around, the systematic overestimate of the final spin is of the 
order of the signal. This systematic overestimate is removed 
if the TTT growth is stopped at about 56-70 per cent of the 
turn-around time, or about 23-35 per cent of the collapse 
time, depending on how these times are defined. However, 
the random scatter from halo to halo is of the order of the 
signal itself. Part of this scatter is due to real non-linear 
evolution, reflected in an uncertainty in the choice of the 
effective time at which the TTT growth should be stopped, 
and part is due to the approximations involved in applying 
TTT already at the initial conditions. The standard trun- 
cation of the potential and smoothing of the shear tensor is 
close to the best one can do in TTT. 



Figure 7. The distribution of angles between model and true 
spin directions in the linear regime, 2 = 50. The three histograms 
(slightly shifted horizontally to improve readability) correspond 
to the different methods of computing the shear tensor (§^). The 
dashed line is the expected distribution for random orientations. 
The misalignment is due to the truncation of the Taylor expansion 
of the potential. Additional error is added by the smoothing on 
the proto-halo scale. 
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Figure 8. The distribution of the angles between model and 
true spin directions in the quasi-linear regime, at 2; = 3. The thick 
histogram is for the standard, truncated and smoothed TTT. The 
thin histogram (which has been slightly displaced to the right) is 
for the angular momentum derived from the Zel'dovich velocities, 
i.e., its predicted direction at any time is the true spin direction 
at 2 = 50. The deviation of this Zel'dovich histogram from perfect 
alignment reflects weak non-linear temporal variations in the spin 
directions. The dashed curve is the distribution of angles in the 
case of random orientations. 
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Figure 9. Evolution of true spin direction. Shown are tlie aver- 
age and the 68.3 per cent confidence interval over the haloes of the 
angle between the angular momentum at time t and its counter- 
part evaluated either at the present time {tg , top panel) or at the 
initial conditions (iso, bottom panel). Linear TTT predicts per- 
fect alignment. The dashed lines show empirical functional fits: 
41°(1 - t/to)°-''^ and 41° {t/to)°-''^ respectively. 



6 SPIN DIRECTION 

The fact that the TTT prediction of the spin amplitude is 
only good to order of magnitude raises doubts about the 
prospects of TTT as a predictor for the direction of the 
spin vector, because a random error of a few tens of per- 
cents in each component is likely to be associated with a 
non-negligible loss of information about the direction. We 
compare the predicted and true directions of the spin vec- 
tor for each halo via the distribution of the angle 9 between 
them over the haloes, proceeding from the linear regime at 
z — 50, to the quasi-linear epoch z = 3, and to the present 
non-linear epoch z = 0, as we did for the spin amplitude. 

In Fig. we compare the three methods of comput- 
ing the shear tensor (94) at the initial conditions, z = 50. 
Zel'dovich velocities would have corresponded to a Dirac S 
function at 9 — 0° , while a random distribution of relative 
directions would correspond to the distribution shown by a 
dashed line. The average angle of misalignment in the lin- 
ear regime is 37°, 26°, and 25° for methods 1, 2, and 3, 
respectively. The corresponding scatter (defined, as before, 
as half the 68.3 per cent confidence interval) is ~ 26°, 18°, 
and ~ 18° in the three cases. Other percentiles of these dis- 
tributions are given in Table ^. As in the case of the ampli- 
tude, the sources of misalignment in the linear regime are the 
second-order truncation of the potential and the smoothing. 
Considering only the haloes with 1000 particles or more, we 
obtain only slightly smaller misalignments, with mean val- 
ues of 33°, 24°, 23° and scatter of 24°, 17°, 15° for the three 
methods, respectively. 

We find that considering the third-order term in the 
expansion of the gravitational potential within the proto- 
halo region (method 1, not shown in the figure) results in a 



Figure 10. As in Fig. g, but at 2 = 0. The distribution of angles 
between the true spin directions at ^ = 3 and z = (not shown) 
practically coincides with the thin histogram, which measures the 
non-linear temporal variations in spin direction. 



very small improvement: the average is reduced from 37° to 
34° , and the scatter from 26° to 24° . Considering the haloes 
for which standard TTT does a poor job and mispredicts the 
spin direction by more than 60°, we find that by adding the 
third-order term one improves the prediction in 70 per cent 
of the cases. For the remaining haloes, where TTT is doing a 
relatively reasonable job, the third-order term does not lead 
to a general improvement and it improves the prediction in 
only 54 per cent of the cases. 

It is interesting to check, already in the linear regime, 
whether a good prediction for the spin direction is associ- 
ated with a good prediction for the spin amplitude. For this 
purpose, we split the halo sample into a "good" group and 
a "bad" group according to whether the spin direction pre- 
dicted by standard TTT and the true spin direction differ 
by less or more than 30°, containing 54 and 46 per cent 
of the haloes respectively. For the "good" group, the aver- 
age and 68.3 per cent confidence interval of I/modci/itruc are 
0.91 and (0.58, 1.18), and for the 'bad' group they are 1.00 
and (0.44, 1.44). We learn that erroneous predictions for the 
spin direction correspond to a large scatter in the ampli- 
tude, while more accurate predictions for the direction tend 
to systematically underestimate the amplitude. We find that 
for 31 per cent of our haloes, the TTT estimates lie both 
within 30° of the true direction and differ from the correct 
amplitude by less than 30 per cent. 

The goodness of TTT predictions correlates with the 
halo mean specific angular momentum. For instance, consid- 
ering the 10 per cent of the haloes with the lowest specific 
angular momentum, we find that standard TTT overesti- 
mates the spin amplitude by an average (median) factor of 
1.6 (1.22) with a scatter of ±0.92, and mispredicts the spin 
direction by 67° (60°) ±40°. On the other hand, for the 
10 per cent of the haloes with the highest specific angular 
momentum, we find that the mean (median) and scatter 
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Table 2. Characteristics of the distributions of the angle between the vectors of the spins 
as predicted by the model and as measured in the simulations. Listed are the mean and sev- 
eral percentiles, as well as the scatter as measured by half the central 68.3 per cent region. 
TTT=tidal-torque theory, NB=Af-body simulations, ZA=Zerdovich approximation. 
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are 0.87 (0.83) ±0.24 for the amplitude ratio and 18° (13°) 
±11° for the angle. These statistics vary continuously as a 
function of the halo specific angular momentum. A simi- 
lar behaviour is found using methods 2 and 3 of but in 
these cases TTT does not introduce a systematic error for 
the haloes with the highest specific angular momentum. This 
behaviour seems to suggest that TTT introduces a typical 
angular-momentum error per particle, and good predictions 
can be only obtained when the halo mean specific angular 
momentum is larger than this typical uncertainty. 

In Fig. ^ we test how well the spin maintains the con- 
stant direction implied by TTT during the quasi-linear evo- 
lution, until z — 3. The spin derived from the Zel'dovich 
velocities at z = 3, which is of course parallel to the spin 
at the initial conditions, has an average misalignment with 
the true spin of only 13° with a comparable scatter. This 
small misalignment implies that the spin tends to keep a 
pretty fixed direction during quasi-linear evolution. When 
this small temporal variation of the angle is convolved with 
the truncation and smoothing involved in applying TTT to 
the initial conditions, the misalignment becomes more sig- 
nificantly larger, with an average of 40° and a scatter of 
~ 30° . The mean evolution and scatter of the true halo spin 
direction all the way to z = is summarized in Fig. |^, either 
relative to the direction at 2: = or relative to the direction 
at z = 50. Empirical functional fits are 41°(1 — t/io)'''^^ and 
41° (t/to)"'^^ respectively. We see a significant true misalign- 
ment developing due to non-linear effects. 

Fig. ^ addresses the full non-linear variation of spin 
direction down to z — 0. When the model spin is computed 
from the particles following the Zel'dovich approximation, 
we obtain the true total misalignments due to non-linear 
effects, as in Fig. ^, with an average and scatter of 41° and 
27°. This true misalignment occurs mostly between z = 3 
and 2 = 0: this part of the misalignment (not shown in the 
figure) has an average and scatter of 37° and 24° , almost as 
large as the total true misalignment. 

Finally, the full TTT prediction, with the additional 



truncation and smoothing, yields a mean misalignment with 
the actual spin at z = of 52° with a scatter of ~ 35°. We 
see that most of the misalignment is due to true non-linear 
evolution of spin direction. This implies that no matter how 
exactly TTT is applied, it is doomed to be a limited predic- 
tor of spin direction in the non-linear regime. 

The bottom line is that non-linear evolution causes sig- 
nificant variations in spin direction, which limit the quality 
of linear TTT as an accurate predictor of spin direction. 
Still, there is a significant correlation between the direction 
predicted by TTT and the final spin direction. 



7 SPIN-SPIN CORRELATIONS 

TTT has been used to predict spatial spin correlations as in- 
dicators of intrinsic spatial alignments of galaxies in the con- 
text of interpreting weak gravitational lensing signals (Cate- 
lan et al. 2001; Crittenden et al. 2001). Our results of the 
previous section already indicate that TTT is a poor predic- 
tor of galaxy spin direction. In order to address directly the 
usefulness of TTT as a predictor of spatial spin correlations, 
we compute in this section two-point spin statistics. 

A summary of the relevance to lensing studies is as fol- 
lows. A shear signal, interpreted as a result of gravitational 
lensing due to the large-scale distribution of matter, has 
been detected by several groups (Bacon et al. 2000; Kaiser 
et al. 2000; van Waerbeke et al. 2000, 2001; Wittman et al. 
2000). Lensing shear maps can be used to directly recon- 
struct the mass distribution, free of galaxy biasing. How- 
ever, a significant uncertainty arises due to possible intrinsic 
coherent alignment of galaxy shapes. Indeed, Brown et al. 
(2000) (see also Pen et al. 2000) claim to have detected a 
non-vanishing correlation between galaxy ellipticities in the 
local universe, over a range of angular separations between 
1 and 100 arc-minutes. Since the galaxy samples used for 
these studies are relatively nearby, a lensing interpretation 
is highly improbable, so this signal is likely to reflect intrinsic 
galaxy alignments. 
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The simplest theoretical approach is to link galaxy 
shapes to halo spins by assuming that galaxy discs are per- 
pendicular to the spin vector of their host haloes. The spins 
are expected to be correlated because the tidal field, that 
plays a key role in generating the angular momentum, is ex- 
pected to have a non-negligible coherency over large scales 
(Catelan & Porciani 2001). Alternative models assume that 
the light profile of elliptical galaxies approximate the shape 
of their hosting haloes (Croft & Metzler 2000; Heavens et 
al. 2000). 

We introduce two statistics for measuring two-point 
spin correlations. The two-point correlation function of spin 
directions is simply, 

77(r) = (L(x).L(x + r)> . (15) 

The function 

,72(r) = ([L(x).L(x + r)]^>-i (16) 

measures the correlation of spin axes independent of the 
sense of rotation, and is therefore more directly relevant to 
the intrinsic correlations of shape orientations that mimic 
lensing signal. 

In Fig. ^ we compare the functions ri(r) and ri2{r), 
at different epochs of the simulation, with the prediction 
of TTT as evaluated from the initial density and velocity 
fields.]^ We see in the simulation that at high redshift there 
is a positive spin-spin correlation at separations of a few 
comoving megaparsecs, but the correlation signal on these 
scales is much weaker at z = 0. This is mainly due to non- 
linear effects. On the other hand, the TTT approximation 
maintains the same correlation function at all times, because 
it keeps the spin directions fixed, by construction. At z = 50, 
TTT matches well the simulation data, and since it always 
builds upon the initial conditions, it maintains the same 
correlation signal at late times as well. The linear TTT ap- 
proximation naturally fails to include the non-linear effects 
that weaken the correlations on scales > l/i~^Mpc, and it 
thus overestimates the correlation at late times.^j 

The weakening of the correlation signal on the scales 
shown can be attributed to two non-linear effects. First, the 
members of galaxy pairs tend to get closer due to the evolv- 
ing clustering. Second, the spin directions themselves evolve 
away from the TTT predictions. As a result of the first ef- 
fect, the linear correlation signal simply shifts to smaller 
scales. The interpretation of the second could be that non- 
linear dynamics erases the linear spin-spin correlations on 
all scales, while non- linear halo-halo interactions build up 

§ The correlation functions presented by Heavens et al. (2000) 
and by Croft & Metzler (2000) seem to be less noisy than ours, 
but this is an illusion caused by the fact that they allowed haloes 
smaller than we do, down to 10 particles, while the resolution of 
their simulations is similar to ours. By including more haloes in 
the sample, one may be tempted to assume that the statistical 
errors are reduced, but the error budget must include the very 
large errors associated with measuring spin and shape for each 
halo based on such a small number of particles (see Bullock et al. 
2001). 

^ We also find that adding the third term in the Taylor expan- 
sion of the potential — see Eq. ( [ll[ ) — leaves the prediction for the 
correlation functions almost unaffected (not shown in the figure). 



a new correlation signal on small scales (note, for example, 
that for r ^ 2/i~^Mpc, t] is positive at z = 50, while it is 
not positive at z — 0). We distinguish between these effects 
by also showing in Fig. |ll] the spin-spin correlation func- 
tions as computed using the TTT predictions for the halo 
spins but located at the actual halo positions at z; = 0. Even 
though this model does a better job than linear TTT, it still 
over-predicts the correlations measured in the simulations 
at z = 0. This demonstrates that non-linear halo-halo inter- 
actions and angular-momentum exchanges with additional 
infalling material play an important role in determining spa- 
tial spin correlations on small scales. Using a simulation of 
higher resolution, one can indeed detect significant spin cor- 
relations at 2; = on scales below 0.4/i~^Mpc (Dekel et al. 
2000). 

Since present-day galactic discs have been probably as- 
sembled sometime around z — 1, we also comput e sp in- 
correlations for haloes selected at that epoch. In Fig. [12, we 
show the results for ~ 5400 haloes containing 100 particles 
or more (left panel), and for ~ 68000 haloes containing at 
least 10 particles (right panel). Linear correlations are sig- 
nificantly weakened by non- linear effects already at 2 = 1. 

It is worth mentioning that the amplitude of the cor- 
relation function ?72(r) for haloes selected at < z < 1 in 
our simulation is in good agreement with the observed align- 
ment of galaxy shapes in the local galaxy population (e.g. 
see Figure 1 in Pen et al. 2000). This indicates that spin 
correlations between dark matter haloes might be related to 
the correlations of galaxy shapes. 

In summary: our results show that linear TTT over- 
predicts spin-correlations for non-linear dark matter haloes 
by a factor of a few. Therefore, caution should be applied 
when using standard TTT to predict the effects of intrinsic 
alignments in weak lensing studies. 



8 ANOTHER MODEL FOR SPIN-SPIN 
CORRELATIONS 

An alternative model for spin-spin correlations have been 
proposed by Lee & Pen (2000, LP) and Pen et al. (2001). 
It is also based on linear perturbation theory, but it tries 
to bypass the difficulties associated with determining the 
inertia tensor of protohaloes by assuming that the spin di- 
rections are statistically correlated with the local tidal field 
in a specific way. The basic Ansatz is 

{LiLj\fij) — ^ Sij — afikfkj , (17) 

where fij — Tij /[TikTki]^^^ is the trace-free unit shear ten- 
sor, and a is a correlation parameter, whose value is to be 
determined. The average is over all possible inertia tensors 
that are compatible with a given realization of the tidal field. 

One then assumes that the inertia tensors of protohaloes 
at different locations are statistically independent (which 
is questionable, since the inertia tensors are found to be 
strongly correlated with the tidal field; see Paper II). This 
allows one to perform the ensemble average in two steps, 
first averaging over the inertia tensors for a given tidal field, 
and then over the realizations of Tij. This leads to 

m{r) oc (T,fc(x)Tfc,(x)f«(x')Ti,(x')), |x' - x| = r. (18) 
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While Tik is non-Gaussian, Tik is assumed to be Gaussian, 
so one can proceed with the calculation if one is willing to 
assume further that 



(f',fc(x)f'fe,.(x)f«(x')f,,(x')) 



(19) 

(r,fc(x)rfc,(x)r.i(x')r,,(x')) 



Applying Wick's theorem, LP obtained this way the approx- 
imation 



V2{r) 'r^ — 



(20) 



where ^M{r) and af j — ^m(O) are respectively the two-point 
correlation function and the variance of the mass overdensity 
field smoothed on the typical proto-halo scale. By including 
all the terms of the same order as (,%f, Catelan & Porciani 
(2001; equation 20) provided the more complete expression: 



9a^ 
4< 



8 



14 



5 9 



14 2 

4 J3.mJ5,M + —J5 M 

5 



(21) 



where Jn,M{r) = nr " dqq" Ca/('?)- For scale-free 
power spectra, 772 is indeed proportional to ^1/, as in 
Eq. (§). 

Spin-spin correlations as predicted by Eq. (Ell) are com- 
pared to the simulation results in Fig. |ll| and Fig. |l^. We 
see that there is a reasonable qualitative agreement, but 
not in quantitative detail. The smoothing length adopted 
to compute the analytic prediction corresponds to the mass 
of the smallest (and most abundant) haloes in each figure. 
Note that, because of the af j term at the denominator, the 
analytic result for 772 (with a fixed a) depends strongly on 
M for separations much larger than the smoothing length. 
The value of the parameter a is fixed by matching the pre- 
diction and the numerical result at r ~ 3/i~^Mpc. For the 
haloes selected at z = 0, it is a = 0.56to o4 at z = 50, and 
a = 0.26'''q q| at z — 0. However, a direct estimation of a 
in Eq. ([l^), using the L — T correlations found in the same 
simulation in Paper II, yields different results. For example, 
considering haloes selected at 2: = 0, we find a — 0.28 ± 0.01 
at z = 50, and a = 0.07±0.01 at 2 = 0, namely, a difference 
by a factor of two to four. This difference seems to indicate 
that either Eq. (|l7|) or some of the simplifying assumptions 
leading to Eq. (IfW are not good approximations. 

It is also worth noticing that, for a g iven spin correla- 
tion function, the value of a in Eq. (Ell) must depend on 
the scale over which the tidal field is smoothed. We find 
that this dependence is likely to be spurious, because when 
we evaluate a directly from Eq. ( p^ ) using two alternative 
smoothing kernels (either corresponding to the halo mass or 
8 times bigger) we see no significant scale dependence (Paper 
II, §7). The dependence on the smoothing radius implied by 
Eq. (pil) probably originates from the questionable assump- 
tion made in Eq. (|l9[). This ambiguity regarding the scale 
dependence may affect the interpretation of the observed 
tidal field, which can be measured reliably only on scales 
much larger than the galactic scales where a is defined. For 
example, Lee & Pen (2001) measure a ~ a few x lO"'^ from 
the PSCz data smoothed by Wiener Filter on large scales, 
and then obtain a ~ 0.2 by extrapolation to galactic scales 
using the ambiguous scale dependence implied by Eq. (|20|). 



By construction, Eq. (^l|) does not account for non- 
linear effects, which affect both the spin direction and the 
clustering growth, but we see that its accuracy is limited 
even in the linear regime. The fact that Eq. (^) with 
a ~ 0.1—0.2 provides an approximate description of spin cor- 
relations in our simulation at < 2 < 1 and in observational 
data should therefore be attributed to a lucky coincidence. 



9 CONCLUSION 

We have evaluated the performance of linear tidal-torque 
theory in predicting the spin of galactic haloes using ~ 7300 
well-resolved virialized dark-matter clumps extracted from 
a cosmological A'^-body simulation. We defined haloes in to- 
day's density field using the standard friends-of-friends al- 
gorithm, but it would be useful for future work to check 
robustness to alternative halo finders as well as to different 
cosmological scenarios. 

We found that, for a given proto-halo at the initial con- 
ditions, TTT provides a successful order-of-magnitude esti- 
mate of the final halo spin amplitude. The TTT prediction 
matches on average the spin amplitude of today's virialized 
haloes if linear TTT growth is assumed until about to/S, or 
about 5/9 of the turn-around time for each halo. The ran- 
dom error, from halo to halo, is about a factor of two. This 
makes TTT useful for studying certain aspects of galaxy 
formation, such as the origin of a universal spin profile in 
haloes (Bullock et al. 2001b; Dekel et al. 2001), but only at 
the level of average properties. 

Non-linear evolution causes significant variations in spin 
direction, which limit the accuracy of the TTT predictions 
to a mean error of ~ 50°. Furthermore, spatial correlations 
of spins on scales > l/i~^Mpc are strongly weakened by 
non-linear effects. This limits the usefulness of TTT in pre- 
dicting intrinsic galaxy alignments in the context of weak 
gravitational lensing (Catelan et al. 2001; Crittenden et al. 
2001). This situation may improve if the orientations of to- 
day's discs were determined by the halo spins at a very high 
redshift, which are better described by TTT. On the other 
hand, we know this only for the haloes selected at 2; = 0, and 
not necessarily for the subhaloes which host disc formation 
at higher redshifts. A detailed study of this effect would re- 
quire a high-resolution simulation in which a detailed galaxy 
formation scheme is incorporated. A preliminary analysis 
performed by using dark matter haloes selected at 2: = 1 
showed that, also in this case, non-linear effects decrease 
linear spin-spin correlations by a factor of a few. 

For practitioners of TTT, we found that the standard 
approximations made in TTT, such as the second order ex- 
pansion of the Zel'dovich potential, and the smoothing of 
the shear tensor on the proto-halo scale, are hard to im- 
prove upon. The inaccuracies in the TTT predictions are 
dominated by real non-linear effects rather than the above 
approximations. 

In order to deepen our understanding of how TTT actu- 
ally works, we investigate in paper II the cross-talk between 
the proto-halo inertia tensor and the external shear tensor. 
We find to our surprise that they are strongly correlated, 
in the sense that their minor, major and medium principal 
axes tend to be aligned, in this order. This means that the 
angular momentum, which plays such a crucial role in the 
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Figure 11. Spin— spin correlation functions, as defined in Eq. (gj) (left) and Eq. (h^) (right). The time-independent TTT predictions 
(open squares) are compared with the true simulation data (filled circles) at 2 = 50 (top) and 2 = (bottom). The effect of allowing for 
clustering evolution (CE) by associating the TTT predictions with the halo positions at 2 = is also shown (stars). The TTT predictions 
were artificially displaced by ±0.03 in log(r) for clarity. Error-bars denote one-standard-deviation statistical uncertainties. The spatial 
separation r is in comoving units 



formation of disc galaxies, is only a residual which arises 
from the little, ~ 10 per cent deviations from perfect align- 
ment of T and /. We find there that the T ~ I correlation 
induces a weak tendency of the proto-halo spin to be per- 
pendicular to the major axis of T, but non-linear changes in 
spin direction erase almost any memory of the initial shear 
tensor, and therefore observed spin directions cannot serve 
as very useful indicators for the initial shear tensor (cf. Lee 
& Pen 2000). 
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On the other hand, the strong T — I correlation investi- 
gated in Paper II provides a promising hint for how to solve 
a long-standing problem in galaxy formation theory, of iden- 
tifying the boundaries of proto-haloes in cosmological initial 
conditions (Porciani, Dekel & Hoffman, in preparation). 

In our studies of the cross-talk between the different 
components of TTT, we have also realized another surprise 
that leads to a revision in the standard scaling relation of 
TTT (While 1984). We find in Paper III (Porciani & Dekel in 
preparation) that the off-diagonal, tidal terms of the defor- 
mation tensor, which drive the torque, are weakly correlated 
with the diagonal terms, which determine the overdensity 
at the proto-halo centre. The latter enters the TTT scaling 
relation via the expected collapse time of the proto-halo, 
and the weak correlation with the torque leads to a modi- 
fication in the scaling relation. The revised scaling relation 
can be applied shell by shell, together with an extended 
Press-Schechter recipe for merger history, in order to ex- 
plain the origin of the universal angular-momentum profile 
of haloes (Bullock et al. 2001b; Dekel et al. 2001). It can also 
be very useful in incorporating spin in semi-analytic models 
of galaxy formation (Mailer, Dekel & Somerville 2002). 
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